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Abstract 

We consider the D4 — D8 — D8 brane system which serves as ultraviolet completion of the 
Nambu-Jona-Lasinio model, where the only degrees of freedom carrying baryon charge are 
fermions. By turning on chemical potential for this charge one may expect the formation of 
the Fermi liquid ground state. At strong coupling we use the dual holographic description 
to investigate the responses of the system to small perturbations. In the chirally symmetric 
phase we find that the density dependent part of the heat capacity vanishes linearly with 
temperature. We also observe a zero sound excitation in the coUisionless regime, whose 
speed is equal to that of normal sound in the hydrodynamic regime. Both the linear 
dependence of the heat capacity and the existence of zero sound are properties of the 
Fermi liquid ground state. We also compute the two-point function of the currents at 
vanishing frequency but do not find any singularities at finite values of the momentum. 
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1. Introduction 

Understanding QCD at finite baryon density or finite chemical potential is an impor- 
tant and difficult problem. Lattice computations are complicated by the sign problem, 
and one largely has to rely on phenomenological models. At large values of the chemical 
potential asymptotic freedom ensures that the physics is determined by the dynamics of 
the quarks near the Fermi surface. The physics simplifies in the planar limit, since some of 
the perturbative instabilities are suppressed (see [1] and references therein for the recent 
work in this direction.). 

It is interesting to compare the predictions of perturbative QCD with those coming 
from string theoretic models. Consider the holographic model of QCD [2,3] at finite tem- 
perature T and baryon chemical potential fi. An incomplete list of references where this 
setup was analyzed includes [4-15]. The holographic model is by no means equivalent to 
QCD, and in fact lacks asymptotic freedom. In the holographic model it is possible to 
adjust the glueball mass scale to be lower than the meson mass scale. This leads to the 
existence of the deconfined but chirally broken phase. The phase diagram then contains 
a line of first order phase transitions between a chirally symmetric phase at larger values 
of fi and T and a broken phase at smaller values of these parameters [5,6]. In addition, at 
large values of ^ there is a phase with condensed baryons [10,12], which is not what we 
naively expect from perturbative QCD. One may ask whether the physics of the system 
resembles that of perturbative QCD, and, in particular whether one can see signatures of 
Fermi liquid formation. 

The question is actually more general that this^. Generic attractive interactions are 
believed to destabilize the Fermi surface, leading to the formation of a gap and a super- 
conducting ground state. (For a holographic description of superconductivity see [22-31] 
and references therein.) In the 't Hooft limit the perturbative quark-quark interactions are 
suppressed, and there is an opportunity for the Fermi liquid ground state to survive, even 
though the coupling is strong. In this paper we attempt to analyze the fate of the fermion 
matter from the DA — DS strings by studying its responses to small external perturbations. 
Since our main interest is not QCD dynamics and the Yang- Mills degrees of freedom sim- 
ply provide the strong interaction between quarks, we consider the (decompactification) 
limit which corresponds, in the field theoretic regime, to a certain UV completion of the 
Nambu-Jona-Lasinio model [32]. (See [33] for more details.) We turn on chemical potential 

^ Recent work on holographic Fermi systems includes [16-20], see also [21]. 
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which is expected to lead to the formation of charged matter, whose ingredients, at least 
in the field theoretic regime are fundamental fermions. We consider the phase where chiral 
symmetry is restored, which corresponds to the D8 — D8 branes falling into the horizon of 
the black hole. Such a phase is perturbatively stable for a wide range of fi and T. 

We show that the density dependent part of the heat capacity at low temperature 
is linear in T. This is the behavior expected for systems with a Fermi surface, where 
only a fraction of quasiparticles is excited at small temperatures. We also observe that at 
arbitrarily low temperatures there exists a massless excitation (zero sound) whose speed 
is equal to the speed of normal sound in the hydrodynamic regime. The existence of 
zero sound is also a feature of the Fermi liquid and corresponds to the deformation of 
the Fermi surface away from the spherical shape. We then compute the current-current 
two-point function at vanishing frequency a; = and finite spacial momentum q using the 
holographic techniques. Such a two-point function should be sensitive to the finite gap in 
the distribution function at the Fermi momentum q = Qf- We do not find any singular 
features however. As we discuss below, this might be related to a non-generic dispersion 
relation of the quasiparticles near the Fermi surface. 

The rest of the paper is organized as follows. In the next section we review the 
DA — D8 — D8 system at finite chemical potential, including phase structure and thermo- 
dynamics. We compute the speed of normal sound in the hydrodynamic regime and show 
that the density dependent part of the heat capacity vanishes linearly with T at small tem- 
peratures. We then study small fluctuations at small temperature (coUisionless regime) 
and find the massless excitation (zero sound) in section 3. The speed of zero sound is equal 
to the speed of normal sound in the hydrodynamic regime. In section 4 we compute the 
current-current two-point function and show that no visible singularities are present. We 
discuss our results in section 5. Appendix A contains the derivation of the speed of normal 
sound in the D3 — D7 system. Appendix B is devoted to the derivation of the boundary 
conditions in the limit of vanishing u. 

2. Thermodynamics 

2.1. Review 

In this section we will review the system studied in [33] at zero temperature. The 
configuration we are interested in consists of three kinds of branes: DA, D8 and D8 
branes. The Nc DA branes are extended along the a;°, • • •,a;^ directions whereas the Nf 
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D8 and D8 branes are parallel and located a distance L apart from each other along the 
x'^ direction. The only massless excitations of the fundamental strings stretching between 
the D8 and DA branes are spacetime Weyl fermions. In particular, 4 — 8 strings contribute 
left-handed ql fermions localized at the X4 = intersection while 4 — 8 strings provide 
right-handed qr fermions localized at X4 = -j. Studying this system at strong coupling, 
i.e., A 3> L 3> /g where A is the fivc-dimcnsional 't Hooft coupling constant, requires taking 
the near horizon limit of the L>4-branes. Assuming Nf <^ Nc, the back-reaction of the D8 
branes on the geometry can be neglected. In this case, it suffices to explore the physics 
of the fundamental matter through the DBI (Dirac-Born-Infeld) action governing their 
dynamics as they propagate in the near horizon geometry of the Nc D4 branes. 

In the following, we summarize some of the results of [6] concerning the DA — D8 — D8 
system at zero temperature but non-zero chemical potential. This requires turning on an 
electric field on the worldvolume of the D8 branes. The value of the chemical potential is 
then read ofi^ from the asymptotic value of the gauge field on the brane. 

Let us start by considering the near horizon region of D4-branes 
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Here t is time, z = 1,2, 3 are the three spatial dimensions of the worldvolume of the D4- 
brane, U is the radial direction and dVL\ the metric of the unit four-sphere. VL^^ and a;4 
denote the volume and volume form of the unit four-sphere respectively. The parameter 
in (2.1) is defined as 

i?3 = TTg^NJl = ttA (2.2) 

where A is the 't Hooft coupling constant. Note that here and in the rest of the paper we 
set a' = 1. 

D8 branes propagating in the geometry (2.1) wrap R^'^ x S'^ and their embedding 
profile is specified by a single function X4{U). Assuming chiral symmetry is restored 
duX^ = and the induced metric takes the following form 

3 3 

^48 - (l) ' i-dt' + dxf) + f^Y" {dU' + UW,) (2.3) 



Studying the system at finite chemical potential, requires introducing non-trivial flux 
Fqu 7^ along the brane worldvolume. The action density (DBI) reads 

Sd8 = j dUUi^l-iduAoY (2.4) 

where N = ^Q,aR^ = ^(27r)~^ ^VI^ " ^^^^ ^^^^ made the gauge choice Au = and 
rescaled according to 

Aq^2'kAq (2.5) 
The conserved charge associated with Aq is given by 

Ui , =d (2.6) 
VI - {duAof 

It follows that the electric field on the brane worldvolume satisfies 

We then determine the chemical potential from the asymptotic value of as 

1 /"°° 1 

''=2^./„ '^^»=i^^'i* (2-8) 



where 7 is defined by 



^. 2. The speed of sound 



7 = ^^'p]^^^] (2.9) 



Holography relates the brane action to the grand canonical potential S = F — ^up, 
where F the free energy density, ji the chemical potential and p the charge density. The 
precise identification at zero temperature is Sd% = — S. Evaluating the DBI action on the 
solution (2.7) we deduce that the grand canonical potential for this phase is given by 

^ = -SDS = -'\j^ldi (2.10) 

where d can be expressed in terms of ji through (2.8). Note that to arrive at (2.10) 
we renormalized the action by subtracting the contribution from the configuration with 
straight branes and no flux. 



Knowledge of the grand canonical potential allows us to determine several thermody- 
namic quantities of the system. The charge density for instance can be shown to be 

A"" 

p=-^ = 2TTdM (2.11) 

Another interesting property is the speed of sound. This is given by 

where P and e denote the pressure and energy density respectively while the derivative is 
taken at constant volume and particle number. At zero temperature, the pressure is equal 
and opposite to the grand canonical potential P — — S while the energy density e is equal 
to the free energy and is given by 

e = ^ + ,jip=^-Mld^ = \P (2.13) 



Using (2.12) we find that 



2.3. Specific Heat 



u' = l (2.14) 



The specific heat capacity Cy is an important physical property of matter which often 
reveals the nature of quasiparticle excitations. The specific heat of a fermionic liquid for 
instance, exhibits linear behavior at low temperatures. On the other hand, one finds for a 
bosonic gas in 3 + 1 dimensions that Cy ~ T"^. In the following we will see that the heat 
capacity of the D4 — D8 — DS system varies linearly with the temperature in accord with 
the predictions of the theory of Fermi Liquids [34] . 

The DA — D8 — D8 at finite temperature and chemical potential has been studied 
in [6]. Here, we will briefly review the results of [6] necessary for the computation of the 
specific heat. The relevant background geometry is given by 
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= fl^s ( t: 1 -^4 = dCs = ^^^^4 
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with f{U) = 1 — The temperature T is related to the minimum value Ut of U as 
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(2.16) 



The induced metric on the straight D8-branes is then 
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dsls^i-B] {-f{U)df + dx 



+ u'^dnl 



R) ^ '^"'"^ ■ "-'^ ■ \RJ \J{U) 
while the DBI action for this configuration at finite chemical potential reads 



(2.17) 



Sd8 = j dUU'^y/l-{duAor 



(2.18) 



where V3 denotes the volume of R^. The equation of motion following from (2.18) has the 
exact same form as the one at zero temperature. Eq. (2.7) is therefore still valid, however 
the chemical potential fx is now given by 



fi = 
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with F 
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the standard hypergeometric function. 
Evaluating the specific heat, as well as the charge density, requires knowledge of the 
grand canonical potential S of the system. At finite temperature, the latter is related to the 
DBI action through S = —TSds- As usual however, the action evaluated on the solution 
(2.7) is infinite. Holographic renormalization [35], [36] is then required for a consistent 
removal of the divergences. Nonetheless, for the purposes of calculating the charge density 
dependent terms of the grand canonical potential any renormalization scheme will suffice. 
This is because turning on a non-trivial chemical potential does not introduce additional 
divergences into the action. The appropriate counterterms are thus fixed and independent 
of the charge density or chemical potential of the system^ . 

We therefore proceed to renormalize the action by subtracting the contribution of 
straight branes with no flux. The result is 



E=AfV3 / dU 

JUt 



Vu^Td^ 



l-F 



7 3 

To' To' 



(2.20) 



^ This can be mainly attributed to the flat brane embedding considered in this paper. We 
thank M.Taylor for explaining this to us. 
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We are interested in the density dependent term of eq. (2.20) which we subsequently define 
as 

9 7 n 7 ."^ ,^2 -| 

(2.21) 



AS = -AfVs'^U^F 



2' 10' 10' f/| 



It is now straightforward to compute the charge density p from 



P = 



dji 



2TTNd 



The specific heat on the other hand is defined through 



C =T(^^ 



where the entropy S is given by 



/a AS 



V dT 



Using (2.19) along with the identity 



/aAS\ /^AS\ f M 



m _i 

\dT)^\dd)^\dn), 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



we find that the density-dependent part of the specific heat^ behaves at low temperatures 
like 



167r^ 



(2.26) 



3. Small Fluctuations and Zero Sound 

In this section we will compute the massless excitation coupled to the density operator 
in the DA — D8 — D8 system at strong coupling. This requires analyzing the linearized 
equations of motions which follow from the action describing the dynamics of D8 and D8 
branes. The full action consists of a Dirac-Born-Infeld (DBI) and a Chern-Simons (CS) 
term. 

The DBI part of the action can be expressed as 

Sdbi,D8 ~ J J d^x J dUe-^y/-det[g + J^] (3.1) 
^ The density independent part behaves like T®. 
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where $ is the dilaton, Q the induced metric and the gauge field strength. We will 
consider fluctuations which are independent of the coordinates of S^. Given that metric and 
gauge field perturbations decouple, it suflBces to set G = G^^^ and expand the gauge field 
^0 = ^0°'' + ^0) = and the gauge field strength Tqu = F^^ + Fo^r, J^ij = Fij. Here 
i,j = 0, 1, 2, 3 and the superscript (0) denotes the background values of the corresponding 
fields. Moreover, we set Ae = 0,^ 9 E S'^ and choose the gauge Au = 0. 
The DBI part of the action for the fluctuating fields is then given by 



Sdbiji = —Af J d^x J dUg{U) 



i i<j 

(3.2) 

where the functions g{U), fi{U) are defined as 



. (3-3) 

MU) = h{u) = 



IJ3 ■'-V ^ U^U^ + d? 

Not all of these functions are independent from each other. In particular, /2 = /1/3. 
We now turn to the contribution of the CS term to the action 

Scs = i^ f F4AAAJ^AJ^ = i ^^^^^^\ c f ANT NT (3.4) 
3! jd% 3! J 

where in the last equality we integrated the RR four form F4 over the four sphere. Re- 
call that the gauge field components Ai with / = 0, 1,2,3, t/ are rescaled according to 
Ai^2ttAi. Expanding (3.4) to quadratic order in the fields we arrive at 

Scsji = i^dRi ^A/- j d^x J dUg{U)^^ [A.F^s + A^F^^ + A^F^^] (3.5) 
With the help of (3.5) and (3.2) we can finally write the full action for the D8 branes as 



Sdsji = J d^x J dUgiU) 



i i i<j 

-lAdRi 4^ (A1F23 + A2F31 + 



(3.6) 

It is convenient to express the equations of motion in the momentum space represen- 
tation where 

Am{x^,t) = j ^e^'=.-^IM(/c^r) (3.7) 



Choosing kfj, = (—a;, 0,0, q) we arrive at the following set of equations 

du[9{U)h{U){duAo)] - q^R'g(U)h(U)A^ - qujR'g{U)h(U)A\\ = 
du[giU)iduA\\)]+uj''R^giU)f2iU)Ai\ + quR^giU)f2iU)Ao = 

du[g{U){duA^,i)] + g{U)R'[u;^f2{U) - q^fsiU)]!^,^ - 4dRhg{U)i^A^,2 = 

du[giU)iduA^,2)] + g{U)R^[uj^f2{U) - q^h{U)]A^,2 + MRiqg{U)^-^A^,^ = 

(3.8) 

To fix the residual gauge invariance we additionally impose Gauss law: 

qA[+uhA'^ = (3.9) 

Henceforth we will focus on the longitudinal modes of the gauge field. Note that the CS 
part of the action affected the field equations for the transverse components only. 

Computing the quasinormal spectrum requires working with the gauge invariant com- 
bination 

E = qAo+ujA\\ (3.10) 
We therefore use Gauss Law to express the first derivative of Aq in terms of E 

^0 = 2 (3-11) 

q^ - uj^fi 

and combine the equations for the longitudinal modes in one 

E" +(^- + f , ] E' + R'M^^h - q')E = (3.12) 

\9 hj 

Primes indicate differentiation with respect to the variable U. Note that to arrive at (3.12) 
we also used the identity /2 = /1/3. 

It is convenient to make a change of variables 



y = 2\l— U=-^ (3.13) 



to express eq. (3.12) as 



y 9 fiq- i^Vi / /i 



Dots denote differentiation witli respect to y wliile 
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Here Jx is defined as 



10. ,10 



fi 10 / 



rio,,io 



/i y \1+Jl'^y^^ 



2Ri 



27A 



(3.15) 



(3.16) 



In tfie vicinity of tlie liorizon (3.14) reduces to a Bessel-type differential equation 

(3.17) 

witli general solution in terms of Hankel functions 



•• 3 • , 

E + -E + (jj^E = 

y 



E{y) = A 



H['\y) , ^H['\y) 



y 



+ B- 



y 



(3.18) 



Imposing the incoming wave boundary condition at the horizon [37-40] singles out one of 
the solutions 



E{y) = A 



H['\y) 
y 



(3.19) 



where H^^^ denotes the Hankel function of first kind. In the limit of small frequencies, or 
to be precise for uy <^1, (3.19) further reduces to 



E{y) ~ — + Aa; 



1 1 

- (1 -27 + i7r + 21n2) - - hiuy 



(3.20) 



with 7 the Euler number 7 ~ .5772. 

On the other hand, for sufficiently small uj and q the last term in eq. (3.14) can be 
neglected. This yields 

E = (3.21) 



^ . 3 g f 1 



y 9 fi(f - ^'^fi 

Eq. (3.21) is analytically tractable. Its general solution can be exressed in terms of 
hypergeometric functions as follows 



E{y) = Co + Ciy' 



1 



Near the boundary in particular (3.22) reduces to 

E{y) ~ Co + ^Ci (l - 



10' 2' 10' ^ ^ 



(3.22) 



y 



(3.23) 
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Imposing normalizability translates to Cq = 0. Hence, the spectrum of quasinormal modes 
will be obtained as a solution to this equation. 

In the vicinity of the horizon the behavior of (3.22) is 

E{y)c^{Co + hC^) + ^ (3.24) 

with a and h defined as follows 

5 o;^ 5/2 a;2x ^ 2 ^2x 

Our next step would be to match the near horizon solution (3.24) to (3.20). Observe 
however that a logarithmic term present in (3.20) is absent from (3.24). The apparent 
inconsistency is resolved by computing the first order correction to the solution given by 
(3.24). In the vicinity of the horizon this is a pretty simple task. It essentially requires 
solving the following inhomogeneous differential equation 

E + ^E = -u;^(co + bCi + ^) (3.26) 

y V y J 

The solution for large y is now given by 



E{y) = {Co + bCi) + ^ - -aC,iu' \ny + 0{y') (3.27) 

y ^ 

which we can readily compare with (3.20) to arrive at 
Co = A 



2 b 
uj 



Ci = - (3.28) 



a 



i(l-2^ + ,.-21„|) ^ 

We are interested in the quasinormal mode with linear dispersion relation in the regime 
of small u and q. Imposing Cq = implies 

^ (^1 - 27 + iTT - 2 In I) - '^Ji^q^ + /I^w^ = (3.29) 

Neglecting higher order terms in u) and q finally leads to 

a;^ = ulq^ «o = 7 (3-30) 
o 

This is exactly equal to the (usual) sound mode computed in section 2. 

As a final note, let us consider corrections to the dispersion relation (3.30). Making 
the substitution w— >a; = ±^^g + 5ijO and expanding eq. (3.29) to first order in 5a; yields 
the following solution 

^ 2 

= ^^S^^' - + ^(^') (3.31) 

10 2//^ lO^fi^ 

Notice that the imaginary part of the dispersion relation is of the type —iq^ in contrast to 
both the typical Fermi liquid behavior [41] and what was observed in [42] , [43] . 
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4. Is there a sharp Fermi surface? 

In the previous section we observed a sound-like excitation in the regime of vanishing 
temperature. It is known that Fermi hquids possess such zero sound mode in the colli- 
sionless regime. It is associated with the deformation of the Fermi surface away from the 
spherical shape. The natural question is whether we can observe the existence of the Fermi 
surface directly. From the theory of normal Fermi liquids we learn that the jump in the dis- 
tribution function can indeed be observed as a singularity in the retarded current-current 
Green's function. More precisely, the two point function at a; = behaves like 

G(.=0,,)~(^-l)log(^-l) (4.1) 

as q approaches the value of twice the Fermi momentum qp. In the previous section, we 
focused on the massless excitation, but the equations can be easily adapted to the case of 
a; = 0, g finite. 

In this section we analyze equation (3.12) for a; = 0. It will be convenient to use the 
variable x = Jiy, where Jl is defined in (3.16) and introduce the rescaled momentum g via 
q = q/Ji. Then eq. (3.12) takes the form 

dlE + - J d^E "L^E = (4.2) 

In the near-horizon region, a; ^ 1 the solution of eq. (4.2) is a linear combination of the 
two solutions, 

£;(!) ^ 1^ £;(2) ^ 1 (43) 

We would like to argue that it is the second solution which is physical. Indeed, as the value 
of q is taken to zero, fiuctuation of the electric field on the brane corresponds to the 
infinitesimal change of the value of //. Yet, in this static situation the boundary condition 
at the horizon imposes E[x^oo) — 0. In fact, in Appendix B we show that E{x^oo) = 
condition smoothly connects to the incoming wave boundary conditions at the horizon, as 
CO is taken to be nonzero. 

We can now use initial conditions at a^rnaa; — 

to integrate eq. (4.2) numerically all 
the way to the boundary. Near the boundary x = the solution of (4.2) is given by 



E = AFi{x)+BFii{x) 
12 



(4.4) 




2 



3 



4 
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Fig 1. Ratio A/B [see eq. (4-4)] a function of q. 



where 



Fi{x) = 



1 - 



~2 2 



2 




) 



+ ..., Fn{x)=x^ + ... 



(4.5) 



The two-point function is proportional to the ratio A/B, which can be extracted from the 
numerical solution with the help of (4.5). The result is shown in Fig.l. As evident from 
this figure, we do not observe any characteristic structure in the wide region of q. Note 
that this is the phase which exhibited such characteristic features of the Fermi liquid as 
linear heat capacity and zero sound. We discuss the significance of these observations in 
the next section. 

5. Discussion 

In this paper we consider the DA — D8 — D8 system in the limit R^^oo which corre- 
sponds to a certain UV completion of the NJL model at weak coupling. Chemical potential 
for the fermions is turned on, and from the field theoretic point of view formation of a 
Fermi liquid is a realistic possibility. The holographic description is applicable at large 
coupling, and a phase transition from the Fermi liquid to some other ground state may 
in principle happen as coupling is varied. However we observed a couple of features that 
suggest that such a phase transition does not happen and we are dealing with a strongly 
coupled Fermi liquid-type system. 

The most significant indication that this is the case is the linear behavior of the density- 
dependent part of the heat capacity at low temperatures (2.26). This is a very generic 
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feature of the systems which have a step-hke distribution function at zero temperature. At 
small temperature the number of excited states is proportional to T, as well as the average 
energy of the excitation, leading to £^ ~ and ~ T. Note that this behavior is in 
sharp contrast with that of Z)3 — Dl system at small temperature, where Cy ~ [42] . 
This difference may be related to the existence of charged bosons in the case of the D3 — D7 
system. Their condensation is presumably responsible for the different low temperature 
behavior of C„. It is impressive that the physics seems to vary smoothly between the field 
theoretic and holographic regimes. 

In section 3 we found a sound mode at vanishing temperature, the zero sound. Again, 
such a feature is characteristic of Fermi liquids, where zero sound is associated with the 
deformation of the Fermi surface away from the spherical shape. It is interesting that the 
speed of zero sound, given by (3.30) is equal to the speed of normal sound propagating 
due to the fermion matter in the hydrodynamic regime [see eq. (2.14).] Note that the 
speed of this "brane sound" differs from that of the sound propagating on the stack of DA 
branes by a factor of V^. (The general formula for the speed of sound on the Dp branes 
is Vg = {5 —p)/{9 —p), see e.g. [44-46].) This is not a surprise since the fermionic degrees 
of freedom are localized on a 3+1 dimensional defect. 

In fact, as we show in appendix A, the value of the speed of zero sound coincides with 
the one of first sound for the D3 — D7 system as well. There, the two expressions match 
as functions of the ratio between the hypermultiplet mass and the chemical potential. The 
fact that the speed of sound in these two regimes is the same is in principle consistent with 
the Fermi liquid theory. As reviewed in [43], when only the first two Landau parameters 
Fq and Fi are turned on, the ratio of zero to first sound velocities goes to one as Fq and 
Fi are taken to infinity. For a generic large interaction F{'d), the ratio still goes to a 
constant. We leave the study of interpolation of the sound mode between the coUisionless 
and hydrodynamic regimes to future work. 

It is interesting that the damping term in the dispersion relation for the zero sound 
behaves like —iq^. This should be contrasted with the —iq^ damping observed in [42], [43]. 
Possible microscopic explanation of this behavior might be related to the non-linear dis- 
persion relation near the Fermi surface. In fact^ this may also explain the absence of 
the apparent singularity in the current-current two-point function at a; = 0, computed in 
section 4. Indeed, the absence of the singularity, which one would generically expect at 

^ We thank Edward Shuryak for pointing this out. 
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q = 2qF is one of the main puzzles that we encounter. It would be nice to have a phe- 
nomenological theory that would explain the observations of the present paper. It is clear 
that such a description must account for strong coupling, since the relevant dimensionless 
parameter A/x is necessarily large for the DBI analysis to be applicable. One obvious dif- 
ficulty in applying Landau's theory of Fermi liquids is the assumption that the number of 
particles is conserved as the strength of the interaction is varied. It is not obvious that 
this should be the case here. 

Finally, we would like to make a comment about the dimensions of the D-branes. 
In this paper we considered the D4 — D8 — D8 system, which has problems with the 
dilaton tadpole beyond the probe brane approximation for the D8 branes. However we 
expect the physics to be similar for other D-brane systems, for example obtained from the 
D4 — D8 — D8 by the T-duality. It would be interesting to investigate these cases in more 
detail. 
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Appendix A. D3-D7 system and speed of sound. 

In this appendix we compute^ the speed of sound for the system studied in [43]. As 
already mentioned in section 2, the speed of (normal) sound in liquids is defined as 




(A.l) 



with P the pressure and e the energy density. Evaluating the speed of sound is then a trivial 
exercise once the grand canonical potential of the system is known. Within the framework 
of gauge/gravity duality, the grand canonical potential is identified with the DBI action 
evaluated on the specific configuration. Here we are interested in the black hole embedding 

^ After this Appendix had been written we received [47] which contains similar result. 
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of the D3-D7 system at zero temperature. Analytic results in this case, including evaluation 
of the action, are given in [48] . According to [48] , but in the conventions of [43] 

where A/", 7, and m are defined in section 3 of [43]. 

At zero temperature the pressure is equal and opposite to the grand canonical potential 
P — — S while the energy density is equal to the free energy and given by 

e = S + //p = -Sdi + /x^l^ = ^AA7(/I2 - m2)(3/l2 + m^) (A.3) 

Given that pressure and energy density are functions of the chemical potential only^ we 
can rewrite (A.l) as 

(A.4) 



and using (A. 2) as well as (A.3) we conclude that 



~2 ~ 2 

2 /i — m 



= ^ (A.5) 

Observe that the speed of zero sound computed in [43] is exactly equal to (A.5). 



Appendix B. The boundary condition at the horizon in the Umit a;^0. 

In section 4, we numerically computed the retarded Green's function for the density 
operator and investigated its behavior for different values of q. Given that the singularity 
which signals the presence of a Fermi surface sits at a; = 0, the calculation was performed 
by setting a; = in (3.14). However, the appropriate boundary conditions at the horizon 
when a; = are not known. Here, we will show that in the limit a;— >0 the incoming wave 
boundary condition reduces to 

lim E{y) = (B.l) 

Let us investigate the region close to the horizon henceforth defined as the regime 
where Jiy ^ 1. Given that we want to study the behavior of the solution at very small 
frequencies, we will not require ^Ji^^y^^ to be large. Eq. (3.14) in this case reduces to 

E+\- + — ^ ]e + uj^\1-^^ ]e = (B.2) 




Note that we take the hypermultiplet mass m as a fixed parameter of our ensemble. 
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Observe that when > 1, (B.2) further reduces to (3.17). In this case, the 

incoming wave boundary condition singles out one of the solutions, given by (3.19). If in 
addition uy <^1, (3.19) behaves like 



E{y) ~ ^ + Auo'^ 

yl 



^ (1 -27 + m + 21n2) - ]-\nujy 



(B.3) 



2 

On the other hand, when ^fi^^y^^ -C 1, eq. (B.2) becomes 



1 

E+— E = (B.4) 

y 

This is precisely (3.14) in the vicinity of the horizon for a; = 0. Its general solution is 

expressed as 

Eiy) = -^ + F (B.5) 

Clearly, one of the solutions vanishes for large y while the other behaves as a constant. 
Which of the two corresponds to (B.3)? 

It is not difficult to see, that there exists an additional, intermediate region where 
(B.2) is analytically tractable. Consider for instance the behavior of the last term for both 
large and small y. 

» 1 

a;M 1 - , \ (B.6) 



^10^10 q2 



Jl'^'y'"' < 1 



It is obvious that this term is negligible so long as^ 



^ V 4 « 2/ « - (B.7) 



Consequently, eq. (B.2) reduces to 



^+|- + — T \E^{) (B.8) 



^ Here we tacitly assumed that > 1- If this is not the case, eq. (A. 8) extends its regime 

of validity in the region <C 1 as well. Repeating the same analysis we can easily see that 

the final result remains unaltered. 
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with general solution 



For large y we can neglect the term which behaves like y ^'^ to get 



E{y)^^+C (B.IO) 

y 

Similarly, in the region of small y (B.9) becomes 

-5^4-10^ + C (B.ll) 

In principle, we should proceed to match eq. (B.IO) to (B.3) and eq. (B.ll) to (B.5). 
Observe however, that (B.IO) does not contain the logarithmic term present in (B.3). To 
recover the logarithmic behavior we need to solve (B.8) in the near horizon region to next 
order in uj. We should therefore consider the following equation 

E+-E = u''(^^ + C] (B.12) 

y \y J 

Solving (B.12) yields 

E{y) = 4 + C - \bu;^ Iny + 0{y^) (B.13) 

y 2 



and combining all the above we deduce that 

(^1 -27 + i7r-21n|) (B.14) 



D 4 



It is clear, that (B.14) vanishes when a; = and the other parameters are kept fixed. 
Note that taking the limit a;— >0 is completely justified given that the matching technique 
is valid in the regime 



4 



a; < /I ( ^ ) (B.15) 
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